Vortex polarity switching by a spin— polarized current 
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The spin-transfer effect is investigated for the vortex state of a magnetic nanodot. A spin current 
is shown to act similarly to an effective magnetic field perpendicular to the nanodot. Then a 
vortex with magnetization (polarity) parallel to the current polarization is energetically favorable. 
Following a simple energy analysis and using direct spin-lattice simulations, we predict the polarity 
switching of a vortex. For magnetic storage devices, an electric current is more effective to switch 
the polarity of a vortex in a nanodot than the magnetic field. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 72.25.Ba, 85.75.-d 



The control of magnetic nonlinear excitations (domain 
walls and vortices) via an electric current is of special 
interest for applications in spintronics jl], The spin- 
transfer effect, theoretically predicted by Slonczewski || 
and Berger and experimentally verified in || ^, [t]] , is 
very promising for magnetization switching with possible 
applications in magnetic storage devices. A perspective 
candidate for this program is the vortex (ground) state 
of a disk-shaped nanoparticle (nanodot) , which provides 
high density storage and high speed magnetic RAM [|| . 

In this Letter we apply the spin-transfer effect to a 
nanodot in the vortex state. Using a continuum descrip- 
tion for the magnetization dynamics we show that a spin 
current mainly acts like an effective magnetic field, ap- 
plied to the sample along the z-axis. The influence of 
a perpendicular field on the vortex properties was stud- 
ied in Refs. fi for easy-plane magnets. There it 
was shown that the vortex behavior depends on its core 
magnetization (polarity) . A vortex with a magnetization 
pointing in the field direction (light vortex) has less en- 
ergy than a vortex with a magnetization pointing away 
from the field direction (heavy vortex). For a strong 
enough field the heavy vortex loses its stability jlO| and 
switches into a light vortex as confirmed by Monte Carlo 
simulations [ pd| . The qualitative similarity between the 
current torque term and the magnetic field suggests that 
an appropriate current can switch the polarity of a vor- 
tex. Using a simple energy argument we show that the 
vortex polarity can be irreversibly switched by a spin- 
polarized current. These analytical predictions are con- 
firmed by lattice spin dynamics. Since an electric current 
is easier to apply and less energetic than a magnetic field, 
we show that this is a more effective way of flipping the 
polarity of a vortex than the magnetic field. 

To realize the program we use a pillar structure, first 
proposed in Ref. [|l2| (see Fig. |l|). An electron current is 
injected in FMi in z-direction, where it is polarized along 



z. The second ferromagnetic layer FM2 is separated from 
FMi by a narrow (about 2 — 3 nm) nonmagnetic layer 
NM JlJ, [H| , so the spin polarization of the current is 
conserved when it flows into FM2 . A spin torque, acting 
on spins of the free layer FM2, causes their precession 
around z with a constant out-of-plane component 0] . 




vortex 

FIG. 1: (Color online). Schematic of the heterostructure. 

The magnetization dynamics in the thin free layer FM2 
is described by the Landau-Lifshitz-Gilbert equation 
with account of the spin-torque effect by a Slonczewski- 
Berger term ||, This equation in dimensionless form 
is 



dt 



[S n xF n ]-e 



S n x 



dSr, 



dt 



■ 3 a(r]) 



Sn_ X 



Sri, X 



1 + b(r]) S n ■ cr 



(i) 



where S n = i^n^n^n) ^ s a classical spin vector with 
fixed length 1 on the site n of the lattice, e is a damp- 
ing coefficient. The effective magnetic field is given by 
F n = —d'K/dSn, where "K is the magnetic energy of the 
free layer. We consider the model of a classical easy-axis 
ferromagnet described by the following Hamiltonian 
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The summation runs over nearest-neighbor pairs (n, n'), 
< S < 1 is the anisotropy constant, "K^d is the mag- 
netic dipole-dipole interaction. The unit vector <x gives 
the direction of the spin polarization ( along z in our 
case), j = J e /Jp is the normalized spin current, J e is the 
electrical current density, J p — /io-Mj \e\d/h, where d is 
the disk thickness, e is the electron charge, /io is the vac- 
uum permeability, M$ is the saturation magnetization. 
The functions a(r]) and b(rj) have the form [|| 
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where < ?/ < 1 is the degree of spin polarization. 

In weakly anisotropic magnets with <5 <C 1, the charac- 
teristic size of excitations £g = 1/ \Jz~8 is larger than the 
lattice constant a — 1 (z is the number of nearest neigh- 
bors in the lattice), so that one can use the continuum 
approximation for the Hamiltonian (0) 
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where Sq is a constant, Udm is an effective de- 
magnetization field which stems from the dipole- 
dipole interaction, and the parameterization S = 
(sin 9 cos 4>, sin 9 sin <fi, cos 0) is used. It is well known 
fill |l6f that a flat cylindrical (disk-shaped) nanodot with 
a radius R smaller than a critical value Rq is in an in- 
plane mono-domain state, whereas for R > Rq the dot is 
in a vortex (curling) state. The vortex state for R > Rq 
becomes energetically preferable as a result of the compe- 
tition between the exchange and the dipole-dipole inter- 
actions. Typical values are i?o = 50 nm for a permalloy 
nanodisk with d = 15 nm g. Generally, the demagne- 
tization field is a nonlocal functional of S, but for flat 
dots, i.e. for dots with the aspect ratio d/R <C 1, the 
action of and the magnetostatic field can be presented in 
a localized form |l7], Q 
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where l e is the exchange length ]l9[ |; hence the de- 
magnetization field is equivalent to an effective easy- 
plane anisotropy with a typical magnetic length scale 

Thus in the continuum limit the current driven dynam- 
ics of flat nanodots is described by the set of equations 
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where the energy functional S is given by Eq. (||) with 
the demagnetization field given by Eq. (||) and a = ±1 
gives the sign of the spin polarization of the current. 

Let us discuss the stationary solutions of Eqs. (||) with- 
out and with current first for homogeneous magnetization 
and then for a vortex state. In mono-domain nanodots 
for j = the magnetization is uniform and lies in the 
ccj/-plane. Equation (||) shows that a current causes a 
finite out-of-plane magnetization 
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where the in plane component 4> rotates around the pr- 
axis with the frequency lo. Because |too| < 1 this state 
exists for small j , it is linearly stable for \j al 2 ~eb\ < e. 
Together with this state there is always the fixed point 
6>o = (resp. #0 = tt) which is stable for j at 2 — s b > e 
(resp. ja£ 2 — eb > — e). 

In the vortex state the magnetization structure is de- 
termined by the expressions 
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where r and \ are the polar coordinates in the xy-plane. 
The localized function f(r) describes the core of the vor- 
tex. In the core the magnetization is nearly perpendic- 
ular to the disk-plane and the index p = ±1, which is 
termed polarity p0[ , determines which way the out-of- 
plane structure is oriented. In nanodots for which the 
vortex (||) is the ground state a current leads to a sta- 
tionary state where the in-plane spin component rotates 
with the angular frequency uj. The in-plane and out-of- 
plane fields can be written as 



= ±-+X + wt + ^(r), 
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and from (||) they follow the equations 
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where V 2 = r~ 1 d r (rd r ) . 

Thus, under the action of the current, the azimuthal 
angle </> acquires an additional dependence ip on the radial 
coordinate r roughly proportional to ebmo/i 2 . However, 
for weak damping (s <S 1) and for small spin polariza- 
tion (77 1) this additional dependence does not change 
significantly the spatial dependence of the out-of-plane 
component cos#. In this limit, in the rotating frame of 
reference where </> is replaced by </> = </> + to i, the Eqs. (||) 
take on the usual Landau-Lifshitz-Gilbert form without 
spin-torque term. A spin current in this case is equiva- 
lent to an effective magnetic field applied along the hard 



FIG. 2: (Color online). The energy of the vortex (j9|) for three different fields: mo = (left panel), mo = —0.2 (middle panel) 
and mo = —0.5 (right panel). The cutoff parameter a = 0.264ao and I — 2.89ao. The switching occurs for mo ~ —0.47. 



z-axis. The energy of the system in the rotating frame 
reads: 



S = 



' cos 9 dr 
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and the model under consideration is similar to the so- 
called cone-state model |], [H^. It was shown in Refs. 

[h]] that the magnetic field removes the degeneracy 
with respect to the vortex polarity p. There are light 
(heavy) vortices, where the core is magnetized along 
(against) the field. The light vortex has always lower 
energy than the heavy one |)) . The range of anisotropy 
and magnetic field for which vortices are stable was found 
in Ref. flf0| . The above mentioned qualitative similarity 
between the current torque term and the magnetic field 
suggests that by applying an appropriate current one can 
achieve switching between the vortex states with different 
polarities. It is worth reminding that in the continuum 
limit the vortex states with different polarities are sep- 
arated by an infinite barrier. However, in real systems 
with their discrete lattice structure, the barrier is finite 
fjlf and switching can occur (see e.g. Refs. 22, ^3], |24| 



where switching due to the noise or an ac magnetic field 
was investigated). 

To give some insight into the current-induced switch- 
ing mechanism, we propose here a simple continuum pic- 
ture, where the discreteness effects are modeled by the 
cut-off parameter a, which is of the order of the lattice 
constant oq. We introduce the simple two-parameters 
Ansatz 
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where the variational parameter \x characterizes the out- 
of-plane core magnetization while the parameter W gives 
the core width. Using this ansatz, one can rewrite the 
energy (0) in the form § w ird [(1 - m%) \n(L/a) + £] , 



where 

cx p(-?) 
£ = 2 

o 



tlntdt 

F¥)F¥)" raoEl(e) (9) 



v = n-mo, £ = -^2. 



Here Ei (£) is the exponential integral. Without the spin- 
current, there exists a double-well potential, which pro- 
vides equal minima for /x = p = ±1 and W = £\/2, see the 
left panel of Fig. ^. Such minima exist only if the cutoff 
parameter aft is smaller than some critical value which 
corresponds to the transition between out-of-plane and 
in-plane vortices H] . For the square lattice such a tran- 
sition is realized for S = 0.297 §]} or £ w 0.9I8a ; this 
corresponds to a/£ as 0.288 and a i=s 0.264ao. 

If the current is turned on, one of the vortex states, 
namely the vortex polarized along too, becomes energet- 
ically preferable, see the middle panel of Fig. 0. When 
the current exceeds some critical value, the barrier dis- 
appears, and the switching occurs, see the right panel of 
Fig. 0. In terms of this variational Ansatz, the switch- 
ing phenomenon can be considered as the motion of an 
effective mechanical particle in the two-dimensional po- 
tential (||). For small currents the particle is captured 
in one of the potential minima but when the current ex- 
ceeds a threshold value one of the minima disappears and 
the particle moves to the another one. The vortex core 
width W ~ £y/2(p + m )/(p — m ) increases with the 
current for light vortices and decreases for heavy ones, 
which agrees with results for the static field |)|. The 
critical value of too, where switching occurs, depends on 
the value of a cut-off parameter, a/£: for the parameters 
presented on Fig. 0, the critical value is mo ~ —0.47. 

To validate our analytical predictions, we performed 
spin-lattice simulations. We stress that the standard 
micromagnetic simulations are not adequate for this 
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FIG. 3: (Color online). Amplitude of S z (t) at the vortex 
center from numerical simulations on a square lattice with 
circular boundary of diameter L — 200. Other parameters of 
the system: r\ = 0.25, e = 0.01, S = 0.03; time is measured in 
units of h/J. The critical current j m —0.0033. 



problem, because the switching phenomenon has a dis- 
crete nature, while micromagnetic simulations are a 
discretization of originally continuum Landau-Lifshitz 
equations. We have integrated numerically the discrete 
(spin-lattice) Eqs. ([!]) over square lattices of size LxL us- 
ing a 4th-order Runge-Kutta scheme with time step 0.01 
and free boundary conditions. The spin system is defined 
by a circular border with the diameter L. We have fixed 
an exchange easy-plane anisotropy parameter S = 0.03 
instead of the dipolar term, corresponding to £/ao = 2.89, 
and we used the damping parameter e — 0.01. We have 
considered system sizes in the range L/ao € (40, 200). 

In our simulations we observed switching between 
heavy and light vortices as a result of the influence of 
the current. Typical pictures are presented in Fig. [|. 
The critical current j w —0.0033 is in a good agree- 
ment with our continuum model calculations: the critical 
value of mo ~ —0.47 (see Fig. |]), which corresponds to 
j w -0.00332. 

To summarize, in this work the spin-transfer effect was 
investigated for the vortex state nanodot. The switching 
of the vortex core magnetization was studied analytically 
and confirmed numerically by direct spin-lattice simu- 
lations. Let us estimate the critical current for a nan- 
odot. Using typical parameters for the permalloy disks 
11(7? = 0.26, M s = 640 kA/m, e = 0.01, d « 10 ran), 
one can estimate that J cr « 10 10 A/m 2 . The total current 
is about 10 mA. For magnetic storage devices an electric 
current is then more effective to switch the polarity of 
a vortex in a nanodot than the magnetic field which is 
about the O.lvlT. 
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